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1.  INTRODUCTION^) 

A  model  for  the  on  frequency'  field  aligned  scattering  has  been  developed 
on  the  basis  of  the  work  performed  to  date  by  all  the  members  involved  in  the 
modeling  task.  The  model  is  intended  primarily  for  use  in  systems  applications 
and  its  presentation  format  has  been  adopted  accordingly.  This  paper  constitutes 
an  attempt  to  present  the  model  in  a  form  that  would  best  meet  its  intended 
objective  and  it  is  hoped  that  it  reflects  the  consensus  of  the  modeling  group. 

The  formulation  of  the  model  begins  with  a  physical  description  of  the 
scattering  medium  which  specifies  the  dimensions  of  the  disturbance,  wave 
number  spectrum  of  the  density  fluctuations  and  a  yield  model  relating  the 
fluctuation  intensity  to  the  heater  power  density.  The  characterization  of  the 
medium  is  an  important  aspect  of  the  modeling  task  and  we  strive  in  this 
presentation  to  see  that  it  is  consistent  with  all  the  available  observations  to 
the  best  possible  extent.  It  is  followed  by  a  description  of  the  radar  properties 
of  the  disturbed  volume  where  the  concepts  of  aspect  sensitivity  and  surface  of 
specularity  are  introduced  as  well  as  the  underlying  principle  relating  the  bistatic 
and  monostatic  scattering  geometries.  The  mathematical  component  of  the  model 
is  presented  in  two  parts,  each  representing  a  different  degree  of  sophistication 
in  computing  the  total  radar  cross  section.  The  first  part  is  concerned  with  a 
general  case  where  a  bistatic  geometry  and  a  realistic  yield  model  are  considered 
and  it  involves  a  numerical  approach.  The  second  part,  considering  a  monostatic 
geometry  and  an  idealized  yeild  model,  offers  an  analytic  solution.  This  case  is 
of  interest  since  a  general  bistatic  radar  problem  can  be  reduced  to  an  equivalent 
monostatic  case  and  can  be  evaluated  to  a  first  order  without  resorting  to  a  complex 
computer  program.  Finally  a  schematic  flow  diagram  is  given  for  computer 
adaptation  of  the  model. 

2.  PHYSICAL  DESCRIPTION  OF  THE  SCATTERING  MEDIUM  (U) 

2A.  Dimensions  of  the  Disturbance  (U) 

The  physical  model  of  the  disturbance  generated  by  the  Platteville  heater 
is  shown  pictorially  in  figure  1.  The  disturbed  volume  containing  the  small 
electron  density  fluctuations  responsible  for  the  field  aligned  scattering  of  the 
radio  waves  is  a  diffuse  region  centered  above  the  heater  at  the  height  where  the 
heater  frequency  equals  the  local  ionospheric  plasma  frequency.  The  strength 
of  the  irregularities  falls  off  above  and  below  this  altitude  with  a  Gaussian  scale 
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MAX.  PLASMA  FREQ. 
300  KM) 


HEATER 

Figure  1.  Physical  Model  of  the  Volume  Containing  Electron-Density  Irregularities  Responsible  for'On- 
Frequency"Field  Aligned  Scattering 


length  of  7.  5  Km  (15  Km  between  l/e  points).  The  lateral  dimensions  of  the 
volume  are  set  by  the  width  of  the  heater  antenna  beam:  for  Platteville  heater 
this  sets  the  Gaussian  radius  of  the  disturbed  volume  at  50  Km  (100  Km  between 
l/e  points)  for  a  heater  reflection  height  of  240  Km.  The  disturbed  volume  can 
thus  be  visualized  as  a  flat  pancake  -  shaped  region  with  diffuse  boundaries.  At 
the  center  of  the  disturbed  volume  the  rms  intensity  of  the  electron  density 
fluctuations  is  typically  1  to  1 .  5%  when  the  heater  is  operating  at  its  full  power 
of  1.9  Mw.  The  physical  description  of  the  distrubance  presented  here  has  been 
arrived  at  primarily  by  requiring  the  resulting  scattering  model  to  be  consistent 
with  the  radar  observations. 

2 B .  Wave  Number  Spectrum:  ( U) 


(U)  The  wave  number  spectrum  of  the  density  fluctuations  of  the  scattering 
medium  is  rented  to  the  backscatter  coefficient  as  (Booker,  1956)^: 

a  / AnY  P  (2Kt2Km,  2Kn)  (1) 


In  the  case  of  strongly  aspect-sensitive  scatter,  the  three  dimensional  wave 

2 

number  spectrum  can  be  written  in  separable  form  as  (Minkoff,  1973)  : 

P(K  ,  K  ,  K  )  =  P,  (K  ,  K  )  P  (K  )  (2) 

'  x’  y  z'  -L  '  x’  y  II  z 

Where  the  two  components  are  defined  in  terms  of  the  corresponding  spatial 
correlation  function  as 

<Ky>  Ky)  =  ff  Rx  (x>  y)  •  exp  |-i(KxX+Kyy)^  dx  dy 

*  -io  *  (3) 


«K,>  * 


J  (z)  .  exp  (-iK^z)  dz 


The  density  fluctuations  are  assumed  to  be  axially  symmetric  around  the  magnetic 

field  and  consequently  the  transverse  correlation  function  Ri  (x,  y)  is  dependent 

**"  2  2- 

only  on  the  radial  distance;  that  is  Rj_(x,  y)  =  Rj_(P),  where  P  =  (x  +  y  )2. 

This  simplification  and  the  use  of  cylindrical  coordinates  lead  to  the  familiar 
Hankel  transform  relation: 


P,_  (KjJ  =  2tt  J  p  Rx  (P)  Jq  (PK)  dP 
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where 


K  2  =  K  2  +  K  2 

-L  x  y 


When  we  are  interested  only  in  fairly  small  values  of  <p,  the  magnetic  aspect 
angle,  we  may  use  the  following  approximations: 


,2  .  2 

JL  +  m  sa  1  and  n  x 


Hence 


Pj_  UKl,  2Km)  =  F>_  (2K) 


(6a) 


and 


P  (2Kn)  =  Pn  (2 K*)  .  (6b) 

using  a  Gaussian  spatial  correlation  for  the  density  fluctuations  along  the 
magnetic  field,  that  is  R„  (z)  =  exp  (-Z2/2L2) ,  we  obtain 

Fj,  (2K*)  =  (2tt  )  ^  L  exp  (-2K2L2^2)  (7) 


using  the  equations  (6a)  and  (7)  we  now  seek  to  obtain  a  relation  between 

the  total  radar  cross  section  and  (2K)  in  orc^er  to  determine,  the  transverse 

spectrum  from  the  radar  measurements. 


a  t  (K) 


4tt 


<t  (K)  dv 

V 


4t: 


lN 


ILL  (^),x 


(2tt)?  L  (2K) 
2 


(8) 


,  y,  z)  exp 


dx  dy  dz 


Where  a  z/R  substitution  was  made,  R  being  the  range  to  the  scattering  volume 
from  the  radar. 


The  intensity  of  the  density  fluctuations  in  a  c  ordinate  system  (X,  Y,  Z) 
centered  on  the  reflectrix  above  the  heater  with  Z  axis  vertically  up  is  expressed 


Where  (AN/N)^  is  the  intensity  at  the  origin,  W  and  A  are  the  Gaussian  radii 

0  — - 2 

of  the  disturbance  in  hoiizontal  and  vertical  directions  respectively.  (AN/N) 

is  expressed  in  the  computational  system  (x,  y,  z)  by  means  of  a  coordinate 

transformation.  The  computational  system  is  centered  on  the  specular  surface 

with  its  origin  right  above  the  heater  and  the  z  axis  perpendicular  to  the  surface. 

The  transformation  is  effected  by 

X  =  x  cos  #  +  •  z  sin 

Y  =  y 

Z  =  -x  sin  0  +  z  cos  6  +  D 


Where  6  is  the  tilt  of  thj  specular  surface  from  horizontal  and  D  is  the  vertical 
mismatch  of  the  specular  surface  from  the  reflectrix.  Using  the  above  form  for 
( AN/N)  ,  the  integral  in  Equation  (8)  is  solved  to  obtain  (Pendyala  and  Thome,  1972) 


^(K) 

J. 


=  (  4v  5  W2  R  (AN/N)2 


[cos2H  sin2#] 

AL 


3.(2  K) 

K 


(10) 


Since  (K^),  where  =  2K,  must  satisfy  the  normalization: 


R.  (0)  =  1  =  1_ 

2tt 


r  Kt  Pj  (Kt)  dKT 

*4 


an 


it  follows  from  Equation  (10)  that 


(AN/N)' 


4  [cos2  8  +  (W2/A2)sin2#J 


'N 


2  ttW2  R 


JkZ  <tt  (K)  dk  (12) 


The  intensity  of  the  density  fluctuation  is  calculated  by  using  the  radar  observa¬ 
tions  o'  Oj  shown  in  Figure  2  to  evaluate  the  integral  in  Equation  (12).  The 
numerical  evaluation  of  the  integral  does  not  include  contribution  from  radar 
frequencies  below  20  MHz.  We  obtain 

<iN/N>o  RMS  *  '* 


when  the  following  parameters  are  used. 


54M 

W  =  50  Km 

10° 

A  =  7.  5  Km 

R  =  900  Km 

5 


(WSbP)  (l>|,d 


MEASURED  TOTAL  RADAR  CROSS  SECTION 


RADAH  WAVELENGTH  XR  (M) 


TRANSVERSE  DENSITY  WAVE  NUMBER  SPECTRUM 
(TWO  DIMENSIONAL!  Ky  =k^  +  «y) 

WAVE  LENGTH  XT  (M) 


O.OOOm 


“i - 1 - 1 - 1 — ; — i — i — i — i — 

2  5  10 

WAVE  NUMBER  KT 


Figure  Z.  Measurement  of  Total  Radar  Cross  Section  a. id  the  Derived  Trans 
verse  Wavenumber  Spectrum  of  Density  Fluctuations 


The  transverse  spectrum  derived  from  the  radar  observations  by  means  of 
Equation  (10)  is  also  shown  in  Figure  2.  It  is  a  three  segment  spectrum  such 


(Kt)  =  0.4 


-  0.  8  X  104  Kt"8-  3 

=  0.  8  X  10"2  K  "L  1 


(0.  8  <KT  <  3.  3) 
(3.  3  x  Kt  6.  6) 

(6.  6  4  KT  <  18.23) 


I  he  curve  drawn  through  the  hatched  region  of  the  wave  number  spectrum 
represents  a  double-  Gaussian  fit  of  the  form: 

Pj_(Kt)  =  2Trj^T2  [exp  -  X  T2  ]  (l3) 

+  (1  -V)  T  2  [exp  -  1  K  2  T  21  I 

2  1  2  J  j 

Whereof  is  a  constant  and  and  are  two  dominant  transverse  correlation 
distances  of  the  density  fluctuations. 

-Yield  Model  [(AN/N2  dependence  on  heater  power  density] 

One  of  the  important  elements  of  the  scattering  model  is  the  density  fluctua¬ 
tion  intensity  (AN/N).  A  model  as  to  how  this  quantity  depends  on  the  heater 
power  density,  referred  to  as  yield  model,  is  presented  in  this  section.  There 
are  two  ways  to  arrive  at  the  desired  model  on  the  basis  of  the  radar  observa¬ 
tions  relating  the  total  cross  section  <rT  to  the  heater  power  P...  One  involves 
performing  a  numerical  inversion  of  the  observed  radar  data  to  obtain  the 
required  solution.  The  other  method  is  indirect,  in  that  a  model  will  be  assumed 
to  start  with,  and  tested  for  its  validity  by  synthesizing  the  radar  cross  sections 
and  comparing  them  with  the  measurements.  The  second  method  has  been 
adopted  here  since  it  is  relatively  simpler.  It  is  assumed  in  the  analysis  given 
here  that  the  heater  power  affects  only  the  intensity  and  not  the  spectrum  of  the 
density  fluctuations. 

Radar  observations  from  PS  V  show  a  wide  variability  in  the  radar  cross 
section  yeild  as  a  function  of  the  heater  power.  The  results  may  be  classified 
into  two  types:  one  revealing  a  constant  slope  and  the  other  a  variable  slope 
with  the  heater  power.  Bulk  of  the  measurements  falls  into  first  category  and 
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tney  are  shown  in  the  form  of  scatter  plot  in  figure  3.  The  cross  section  has  a 
power  law  dependence  on  the  heater  power  with  the  exponent  ranging  over  0.  3 
to  1  with  a  median  value  of  about  0.  65.  In  the  case  of  SRI  data,  it  was  found 
that  lower  slopes  are  in  general  associated  with  the  conditions  when  the  heater 
is  operating  close  to  the  critical  frequency.  For  these  meas  ;rrmer.ts  ol 
constant  slope,  the  yield  model  relating  (AN/N)2  to  the  heater  power  density 
is  the  same  as  that  relating  the  cross  section  to  the  heater  power.  In  some 
instances,  the  second  type  of  behavior  with  variable  slope  is  noted  although 
the  conditions  under  which  this  happens  are  not  yet  identified.  Influenced  by 
the  observed  behavior,  we  have  started  with  a  simple  two-step  model  with 
density  fluctuations  increasing  linearly  with  power  density  up  to  a  certain  level 
and  0.  33  times  as  rapidly  above  that  level  (Figure  4).  Based  on  the  model,  the 
computation  of  the  normalized  total  cross  section  as  a  function  of  the  heater 
power  proceeds  as  follows: 

Let  the  heater  antenna  pattern  be  Gaussian  so  that  the  power  density 
distribution  on  the  reflectrix  is 


PD  =  Po  exP  [-(R/W)2] 


Where  Pq  is  the  power  density  at  the  center,  R  is  the  radial  distance  and  W  is  the 
Gaussian  radius  of  the  antenna  beam.  Let  Pg  be  the  power  density  level  above 
which  the  density  fluctuation  power  law  dependence  changes  from  1  to  0.  33.  Wher 
the  peak  power  density  Pq  is  greater  than  Pg,  the  radius  of  the  region  where  the 
change  occurs  in  the  power  law  is  given  as 

RC  =  W^»(Po/Ps)  (15) 

The  distribution  of  the  density  fluctuations  or  the  reflectrix,  when  normalized 
to  1  at  the  center  for  maximum  heater  power,  is 


(AN/N)2  =  exp  [-  0.33  (R/W)2  ] 


(R/.R  ) 
—  c 


=  exP  C-  °-  33  (Rc/W)2]exp  [-(R2  -  R2)/W2  ]  (R>k  ) 
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When  the  heater  is  operating  with  on-axis  power  density  less  than  Pg,  then 
R  =  0 

c  _ 

and  (AN/N)^  =  (Po/Pg)  (AN/N)g  exp£  (R/W)2] 

where  (AN/N)2  is  the  density  fluctuation  value  associated  with  the  power  density 
PS‘ 

In  order  to  compute  the  radar  cross  section,  however,  the  density  fluctuations 
had  to  be  calculated  on  the  specular  surface.  Consider  an  elemental  area  dA  on 
the  reflectrix  at  some  radial  distance  R  and  azimuths  as  shown  in  Figure  5. 

The  corresponding  elemental  area  on  the  specular  surface  dA*  is  given  as: 

dA1  =  RdRd  <*>/cos2  0 


and 


2„ 

cos  0 


2 

cos  <t> 


cos 


■  2 

sin  <p 


where 


<p  =  azimuth  angle 

(I  =  inclination  of  the  specular  surface 


The  density  fluctuation  intensity  tapers  off  from  the  reflectrix  in  a  Gaussian 
manner  with  a  vertical  scale  size  denoted  as  A.  Hence  the  intensity  on  dA 
is  (AN/N)2  [exp  -  (d/A^  ]  where  d  is  the  vertical  distance  R  tan  6  separating 
dA  and  dA7.  The  incremental  radar  cross  section  dff  contributed  by  dA  can  be 
written  as: 


d«r  =  C  (AN/N)p  exp  -  (d/A)2  R  ~1  dR  d<£ 


(16) 


COS  0' 


where 


C  -  A  constant 


Integrating  Equation  (1  5),  we  obtain 

2tt 


°T  -  C 


If 


P3 


o  © 


(AN/Nr  exp  -  (d/A)' 


dR  d 4> 


(1 
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The  above  surface  integral  is  numerically  evaluated  by  dividing  the  surface  into 
concentric  rings  as  shown  in  Figure  5.  The  outermost  ring  encloses  all  the  area 
of  interest  with  dimensions  twice  the  Gaussian  diameter  of  the  disturbance.  The 
cross  section  is  calculated  for  different  power  values  relative  to  its  value  when 
the  heater  is  operating  at  the  maximum  power.  Figure  6  shows  a  comparison 
between  the  synthesi  ,ed  (dashed  curve)  and  the  observed  cross  section  dependence 
on  the  heater  power.  Given  the  uncertainties  iu  the  yield  measurements,  the 
simple  model  presented  here  can  be  ct.nsidtj.ed  adequate  to  account  for  the 
power  dependence  of  the  radar  cross  section.  Some  caution  should,  however,  be 
exercised  in  adopting  a  simple  yield  model  because  of  the  wide  variability  seen 
in  the  measurements. 

3.  RADAR  PROPERTIES  OF  THE  DISTURBED  VOLUME  (U) 

3A*  Aspect  scnsititivy  and  the  concept  of  a  surface  of  specularity  (U) 

The  high  aspect  sensitivity  of  the  scatterers  leads  to  an  important  modeling 
concept  -  the  surface  of  specularity.  The  importance  of  the  concept  is,  first  of 
all,  that  it  helps  to  physically  understand  the  radar  results,  and  secondly,  that 
it  reduces  the  numerical  task  of  computing  radar  scattering  crosj  sections  from 
a  3-dimensional  problem  to  a  2 -dimensional  problem.  Consider  first  the  optical 
analog  shown  in  Figure  7.  The  glass  rods  are  aligned  with  the  earth’ s  field 
and  represent  the  small  field -aligned  fluctuations  in  refractive  index  caused  by 
heating  the  ionosphere.  The  light  from  the  candle  flame  is  partially  reflected 
by  the  glass  rods,  and  an  observer  behind  the  candle  sees  an  image  of  the  flame 
at  the  specular  point  on  each  rod.  If  there  are  many  closely  spaced  rods,  the 
images  from  each  will  blend  together  to  form  a  bright  'specularity  surface*” 
on  which  the  specular  points  for  all  rods  lie.  The  only  part  of  the  rod  that 
matters  for  the  observer  is  the  segment  that  coincides  with  the  specularity 
surface.  The  corresponding  picture  for  the  radar  case  is  shown  in  Figure  8. 

The  glass  rods  have  been  replaced  by  the  field-aligned  electron-density 
irregularities  caused  by  the  heater.  In  this  case,  the  height  where  the  heater 
frequency  matches  the  local  plasma  frequency  is  well  above  the  specularity 
surface  for  the  radar  and  consequently  only  the  weak  lower  tips  of  the  irregularities 
are  ‘seen*  by  the  r^dar.  The  received  signal  strength  will  consequently  be  weak. 
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RELATIVE  RADAR  CROSS  SECTIO 


RELATIVE  HEATER  POWER  (dB) 


Figure  6.  Comparison  Between  Experimental  and  Computed  Radar  Cross 
Section  Dependence  on  Heater  Power 
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OPTICAL  ANALOG  : 

>  RODS  ALIGNED  WITH  MAGNETIC  FIELD 


icept  of  a  Surface  of  Specularity  for  Reflections  from  a  Bundle 
Scatterers  (Glass  rodsl 


The  situation  can  be  improved  by  lowering  the  heater  frequency  (and  consequently 
the  height  where  the  heater  frequency  matches  the  plasma  frequency)  till  the 
disturbed  volume  coincides  with  the  surface  of  specularity.  The  concept  is 
qualitatively  the  same  for  bistatic  as  well  as  for  monostatic  paths:  The  procedure 
is  to  compute  the  surface  of  specularity,  given  the  end  points  of  the  path  and  a 
model  of  the  earth  s  field,  and  then  use  the  strength  of  the  scatterers  lying  on 
the  surface  to  compute  the  effective  scattering  cross  section  of  the  disturbed  volume. 
This  converts  the  computational  problem  from  a  volume  integration  to  a  surface 
integration  with  a  corresponding  reduction  in  complexity. 

3 B .  Bistatic  Versus  Monostatic  Scattering  Geometries  ( U) 

(U)  When  the  receiver  and  transmitter  of  a  communications  path  are  separated 
from  each  other  (bistatic  rather  than  monostatic  operation)  the  scattering  proper¬ 
ties  of  the  disturbed  volume  are  modified.  The  principal  modification  is  that  the 
size  of  the  irregularities  that  support  scatter  over  bistatic  paths  is  larger  than 
the  size  that  supports  backscatter  at  the  same  frequency.  The  situation  is  shown 
in  Figure  9.  The  Bragg  condition  for  constructive  interference  says  that  the 
spatial  Fourier  component,  which  counts  in  determining  the  strength  of  the  wave 
scattered  through  an  angled,  is  the  component  of  spatial  wavelength  s  -\/Z  sin  (o/Z), 
where  A  is  the  wavelength  of  the  incident  radio  wave.  When  this  condition 
is  met,  the  waves  scattered  from  adjacent  crests  of  the  spatial  Fourier  components 
will  differ  in  phase  by  exactly  one  wavelength  and  thus  will  add  coherently  in  the 
direction  of  interest.  Figure  9  Illustrates  the  fact  that  the  spatial  wavelength  of 
interest  becomes  longer  as  the  path  changes  from  monostatic  to  bistatic.  Since 
the  strength  of  the  spatial  frequency  spectrum  increases  monotonically  as  the 
wavelength  becomes  longer  (see  Figure  3,  bottom),  the  disturbed  volume  will 
scatter  more  strongly  over  a  bistatic  path  than  over  a  backscatter  path,  all  else 
being  equal. 

(U)  The  bistatic  gain*  predicted  by  the  model  for  a  90°  scattering  angle  is 
shown  in  Figure  10.  The  shaded  area  shows  the  range  in  which  backscatter 
measurements  of  total  radar  cross  section  fall  as  a  function  of  frequency.  The 
solid  backscatter  "  curve  shows  that  the  model  predictions  are  in  agreement 
with  the  observations,  and  the  solid  '  bistatic ''curve  shows  the  increase  in  cross 
section  expected  for  a  90°  scattering  angle.  Gains  of  10  dB  or  more  can  be 
expected,  depending  on  frequency. 
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igure  9-  Illustration  Showing  that  the  wavelength  of  the  spatial  frequency  component  that  controls  the  st 
of  the  signal  scattered  through  an  angle  9  increases  as  the  scatter  geometry  changes  from  mon 
static  to  bistatic.  (Here,  A  is  the  radiowavelength  and  s  is  the  wavelength  of  the  spatial  Fourie 
component  that  produces  constructive  interference  in  the  direction  of  interest)  (U) 


4.  MATHEMATICAL  MODEL  FOR  o  (U) 


4A .  General  Case  (Bistatic  Geometry  and  Realistic  Yield  Model)  ( U) 

The  formulation  of  the  mathematical  model  presented  here  involves 
basically,  the  derivation  of  an  expression  for  the  radar  cross  section  per  unit 
area  on  the  specular  surface.  The  total  radar  cross  section  is  computed 
by  numerical  evaluation  of  the  surface  integral  of  <r^.  The  radar  cross  section  per 
unit  volume  of  the  scattering  medium,  which  must  be  integrated  along  the 
direction  normal  to  the  specular  surface  to  obtain  c^,  is  expressed  by  Booker* 
as: 

%  =  tt2  sin2X  (AN/N)2  P  |k<£2  -^),  K(m2  -  m^,  Kfrn,  -  n^  | 
XN  (18) 

where 


X  =  Angle  between  incident  E  vector  and  the  direction 

of  scattering 

=  Background  plasma  wavelength 
(Xy  m^n^)  =  Direction  cosines  of  incident  wave  vector 

and 

i 

(^2 ,  m2>n2)=  Direction  cosines  of  scattered  wave  vector 


Resolving  the  three  dimensional  wave  number  spectrum  into  parallel  and  trans¬ 
verse  components  as  before  and  assuming  the  spectrum  to  be  cylindrically 


symmetric,  we  can  write 

1_ 

2  \ 

a  =  tr  2 sin2  X  (AN/N)2  P.  Jk 

[UZ-^)Z  +  (Hl2-mi)2]  I 

V  4 

XN  V 

)  K( n2  ~  ni ) 


(19) 


Now  let  us  consider  a  unit  volume  at  point  P  above  the  elemental  area  dA^  of 
the  specular  surface  as  sketched  in  Figure  11.  The  elemental  area  is  located 
at  (X^,  Y^,  Z^)  in  the  coordinate  system  centered  on  the  heater  reflectrix  and 
the  density  fluctuation  at  this  point  is 

(AN/N)2  ^  (AN/N)2r.  exp  |- ( Zi/ A)2  |  (20) 
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Figure  11.  Sketch  of  3istatic  Radar  Geometry  showing  the  coordinate  systems  on 
the  specular  surface  and  the  heater  reflect  rix  used  in  the  model 
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Where  (AN/N)  is  the  density  fluctuation  on  the  reflectrix  at  a 


2  1 


radial  distance  =  (X^  +  Y.  )  2  and  it  is  determined  by  the  yeild  model 

described  earlier.  The  fluctuation  intensity  at  point  P,  that  is  at  (0,  0,  z) 
in  the  system  centered  on  dA^.,  is  given  as: 

v2 


(AN/N)2  =  (AN/N)j  exp  -  (  z  cos7 


J z  cos*y 
\  AT 


(21) 


Where  y  -  (90  -  I)  and  I  is  the  dip  angle, 


On  the  assumption  of  Gaussian  spatial  correlation,  we  get  for  the  parallel 
component  of  the  density  spectrum 


K  (n^  -  n^)>  =  (2tt  )  2  L  ei 


,p-jKV<n2-n/j 


(22) 


using  the  geometry  shown  in  Figure  12,  the  direction  cosines  of  the  incident 

and  scattered  wave  vectors  can  be  related  to  the  aspect  and  scattering  angles 

t  and  6  as : 
s 

(n0  -  n.)^  =  4  sin2  </>  /2  sin2'/' 

cl  s 


(T2  +  (m2  •  mj),  =  4 


sin  /2  cos  'I' 
s 


In  the  case  of  strong  aspect  sensitive  scatter,  the  effective  scattering  is  limited 
to  small  values  of  ip  and  we  can  approximate  the  above  relations  to 

(n  -  n  )2  «  4 z 2 / P?2  Where  R*  =  2R1  K2  (23a) 

6  1  i  r> 


(Rj  +  R2) 


and 


(^2  --(j)2  +  (m2  -  mj)2*!  4  sin2 


(23b) 


using  equations  (21),  (22)  and  (23a,  b),  we  get  for  (19) 


=  (2it  )  2  it  2  sin2X 

(AN/N)f  L 

P±  (2K  sin  <4  /2) 

4 

l 

an 

1  r  2 

„„2 T  2  T 

2 

exp )  - | cos  7 

+  2K  L  , 

z 

I  L2 

r'  2  J 

(24) 
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Figure  12.  Geometry  showing  the  Relation  Between  the  Direction  Cosine 
Scattering  Angles  f  and  .  (U) 
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Integrating  Equation  (24)  with  respect  to  z  and  letting  Y,  X  •  an^  take 
values  appropriate  to  the  centpr  of  dA^  will  result  in 


U  A  • 

Ai  = 


2  2 

tt  sin  X  . 


.  A2  R^  cos2  Y  . 
1  + _ i 


J  ~~  (AN/N)i  (2K  sin  %,/> 2)  (25) 
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si 


•N. 


The  total  cross  section  is  obtained  by  numerical  evaluation  of  the  surface 
integral  of  Thus 


a  T  =  f° AdA  =  S 


N 


^Ai 


dA. 


(26) 


i  =  1 


where  N  is  the  number  of  elements  into  which  the  active  region  of  the  specular 
surface  is  divided.  Figure  13  shows  a  scheme  for  computer  adaptation  of  the 
scattering  model. 

4B.  Special  Case  (Monostatic  Geometry  and  Idealized  Yield  Model) 


The  scattering  model  for  the  total  cross  section  can  be  presented  in  an 
analytic  form  for  the  special  case  of  a  monostatic  geometry  and  an  idealized 
yield  model.  The  solution  is  of  some  interest  since  a  general  bistatic  radar 
problem  can  be  reduced  to  ar  equivalent  monostatic  case  and  can  be  evaluated  to 
a  first  order  without  resorting  to  a  complex  computer  program.  The  approximate 
nature  of  the  solution  is  a  consequence  of  two  simplifying  assumptions  that 
readily  permit  an  analytic  description  to  the  model.  The  assumptions  are  that  ths 
magnetic  field  is  constant  over  the  dimensions  of  the  scattering  volume  and 
that  the  intensity  of  the  density  fluctuations  is  directly  proportional  to  the 
power  density  on  the  heater  reflectrix.  Starting  with  Equation  (1)  and  following 
the  steps  given  by  Pendyala  and  Thome  (1972)  ,  we  obtain  for  a  pulse  radar 
with  a  Gaussian  pulse  length  p  and  beam  width  B: 


(K)  = 


4  2 

2  it  ’\v  A  R 


ffl 


Mf  P,  (2K) 
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1 

£  cos2  6  +  W2 
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(27) 
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Figure  13.  Flow  Chart  for  Computer  Adptation  of  the  Model  (U) 
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where  M^  is  a  factor  arising  due  to  mismatch  between  heater  reflectrix  and 
specular  surface  and  i?  given  as: 

Mf  =  exp  -  j  D^_  (1  -  A)j  »  exp  -  Jd2/A2J  (28) 


A  is  a  function  of  9  ,  A,  W,  p,  B  and  R  and  is  much  less  than  1  for  most  cases  of 
interest  thus  leading  to  the  above  approximation.  In  the  case  of  strong  aspect 
sensitive  si  attering}5^^  1  further  simplifying  the  equation  for  a  .  By  letting 
the  pulse  length  and  the  beam  width  to  exceed  the  dimension  W  of  the  scattering 
volume,  Cp  can  be  made  to  approach  the  total  cross  section  c^.  On  the  basis 
of  the  relationship  between  monostatic  and  bistatic  radar  configurations  discussed 
earlier,  the  total  cross  section  for  a  bistatic  radar  illuminating  the  entire 
disturbed  volume  obtains  the  form: 


"x  (K) 


2tt4  sin2  X  WA  R; 


Mj.  P  (2K  sin  ^s/2) 


[ 


cosZ  9  +  (W/A)2  sin2 


] 


(29) 


Where  <f>s  is  the  scattering  angle  and  the  equivalent  range  R*  =  ZR^/^  +  R2). 
The  above  equation  will  be  useful  for  the  purpose  of  quick  evaluation  of  a  given 
radar  link. 


5.  MODEL  VERIFICATION  (U) 

According  to  the  model  presented  above,  there  will  be  an  enhancement  in 
the  radar  cross  section  to  be  gained  by  going  from  monostatic  to  bistatic 
configuration  when  operating  in  the  region  where  the  spectral  function  Pj_(Kt) 
decreases  with  the  increasing  wavenumber.  A  multi-station  CW  scattering 
experiment  designed  to  verify  the  model  prediction  was  conducted  by  SRI  using 
an  operating  frequency  of  144  MHz  that  falls  in  the  region  of  steep  slope  of  the 
cross  section  -  versus  -  frequency  curve  (Frank,  et  al ,  1973).  4  The  experiment 
involves  a  comparison  of  the  cross  sections  observed  over  a  near-backscatter 
path  with  <£g  =  168°  (Ocotillo-Lancaster)  and  a  forward  scatter  path  with  4>  =  89° 

(Ocotillo-Huntsville).  The  observations  show  a  10  dB  gain  in  the  cross  section 
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for  the  forward  scatter  path  relative  to  the  backscatter  case.  The  scattering 
model  predicts  an  enchancement  of  about  14  dB  on  the  assumption  that  the 
polarization  and  mismatch  effects  are  small.  The  difference  noted  between  the 
model  prediction  and  the  observation  is  considered  to  be  within  the  error  bounds 
of  the  experiment. 

6.  CONCLUDING  REMARKS  (U) 

The  scattering  model  presented  here  deals  only  with  the  so  called  'on- 
frequency  component  of  scatter  from  the  heated  volume.  There  are  other 
components  (e.  g.  ,  ion  lines  and  plasma  lines),  but  these  are  relatively  weak 
and  appear  to  be  of  secondary  importance  for  applications.  The  model  has  been 
developed  on  the  basis  of  information  available  primarily  from  the  field 
experiments  reported  to  date.  In  order  to  place  the  model  on  a  firmer  ground, 
relevant  laboratory  and  theoretical  contributions  also  need  to  be  incorporated 
into  the  model.  Although  the  model  is  based  upon  the  observations  limited  to  the 
conditions  of  ordinary  mode  heating  of  the  daytime  F-region,  the  underlying 
concepts  apply  as  well  to  any  other  situation  where  the  scattering  mechanism 
is  iuentified  to  be  the  same.  Specifically,  the  model  needs  to  be  extended  for 
the  E  region  and  the  nighttime  heating  conditions. 


APPENDIX 

DETERMINATION  OF  SPECULAR  SURFACE  (U) 

(U)  The  specular  surface  is  defined  by  a  region  of  points  where  the  external 
bisector  to  the  angle  between  the  incidence  and  scattering  directions  achieves 
perpendicularity  to  the  direction  of  the  earth’ s  magnetic  field.  The  determination 
of  the  specular  surface  within  the  volume  of  the  scattering  medium  for  a  given 
radar  configuration  is  one  of  the  main  modeling  tasks.  An  analytic  solution  was 
obtained  to  the  problem  for  the  case  of  curved  earth  and  dipole  model  for  the 
earth’s  magnetic  field  but  neglecting  the  ionospheric  refraction  effects.  The 
solution  provides  the  height  of  a  specular  point  for  a  specified  location  on  earth 
and  it  results  from  a  cubic  equation  which  is  formulated  by  equating  the  angle 
between  the  incident  ray  and  the  magnetic  field  to  that  between  the  scattered  ray 
and  the  magnetic  field. 

The  solution  was  adopted  to  a  CDC  6700.  The  computer  output  for  a  specified 
radar  geometry  and  heater  location  is  a  matrix  of  325  specular  heights  calculated 
at  equally  spaced  intervals  within  a  region  of  180  Km  square  centered  over  the 
heater.  Figure  14  shows  a  specimen  display  generated  for  the  case  of  Los  Lunas  - 
Soccorro  radar  path  and  Platteville  heater.  The  program  can  compute  and 
display  the  ranges  from  the  transmitter  and  receiver  to  each  of  the  325  specular 
points  in  space  if  so  desired.  A  more  general  program  has  been  developed  at 
SRI  to  compute  the  specular  surface  which  takes  the  ionospheric  refraction  effects 
also  into  account. 
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LOS  LUNAS  (34.798°N,  106.66°W)  -  SOCORRO  (33.963°N,  106.889°W) 
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igure  14.  A  Sample  Display  of  Specular  Height  Matrix  above  Platteville  Heater  for  Los  Lunas  -  Soccorio 
Radar  path 
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planning  and  executing  the  USAF  exploratory  and  advanced 
development  programs  for  electromagnetic  intelligence 
techniques ,  reliability  and  compatibility  techniques  for 
electronic  systems ,  electromagnetic  transmission  and 
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